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A detection of the level of non-Gaussianity in the CMB data is essential to discriminate 
among inflationary models and also to test alternative primordial scenarios. However, 
(-H ' the extraction of primordial non-Gaussianity is a difficult endeavor since several effects of 

C*^ , non- primordial nature can produce non-Gaussianity. On the other hand, different statis- 

I 1 tical tools can in principle provide information about distinct forms of non-Gaussianity. 

O ' Thus, any single statistical estimator cannot be sensitive to all possible forms of non- 

}~[ ' Gaussianity. In this context, to shed some light in the potential sources of deviation 

1^ ^ from Gaussianity in CMB data it is important to use different statistical indicators. In 

C^ ■ a recent paper we proposed two new large-angle non-Gaussianity indicators which pro- 

vide measures of the departure from Gaussianity on large angular scales. We used these 
indicators to carry out analyses of non-Gaussianity of the bands and of the foreground- 
reduced WMAP maps with and without the KQ75 mask. Here we briefly review the 
formulation of the non-Gaussianity indicators, and discuss the analyses made by using 
our indicators. 



1. Introduction 
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lO ' A detection of the level of primordial non-Gaussianity in the CMB data is crucial 

^^ , to discriminate inflationary models and also to test alternative scenarios for the 

physics of the early universe. Clearly the study of detectable non-Gaussianities in 
the WMAP data ought to take into account that they may have non-cosmological 
origins as, for example, unsubtracted foreground contamination, unconsidered point 
^\f ' sources emission and systematic errors. Deviation from Gaussianity may also have 

a cosmic topology origin (see, e.g., the review Refs. [Hand related Refs. [21). If, on 
the one hand, different statistical tools can in principle provide information about 
distinct forms of non-Gaussianity, on the other hand one does not expect that a 
single statistical estimator can be sensitive to all possible forms of non-Gaussianity 
in CMB data. In view of this, a great deal of effort has recently gone into verifying 
the existence of non-Gaussianity by employing several statistical estimators (an 
incomplete list of references is given, e.g., in Refs. [SHlJand references therein). 

Recently we have proposecP' two new large-angle non-Gaussianity indicators, 
based on skewness and kurtosis of large-angle patches of CMB maps, which provide 
measures of the departure from Gaussianity on large angular scales. We used these 
indicators to search for the large-angle deviation from Gaussianity in both band 
and foreground- reduced mapa^with and without a KQ75 mask (see also the related 
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Refs. in])- Here we briefly summarize the main results of Refs. |l]and[51 

2. Results and Concluding Remarks 

A constructive way of formulating our non-Gaussianity indicators S and K from 
CMB data is through the following steps: 

(a) Take a finite set of points {j — l,...,A'c} homogeneously distributed on the 
CMB celestial sphere S^ as the centers of spherical caps of a given aperture 7; 
and calculate for each cap j the skewness (Sj) and kurtosis (Kj) by using that 

S^^W^Y.iT^~T^f and K,^j;^-^J2iT^-T,y-3, (1) 



3 i=l ^ 1 i=l 

Vp '" ^ 



where TVp is the number of pixels in the j *'^ cap, Ti is the temperature at the i ^^ 
pixel, Tj is the CMB mean temperature of the j**^ cap, and a is the standard 
deviation. The numbers Sj and Kj obtained in this way for each cap can be seen 
as a measure of non-Gaussianity in the direction of the center [Oj , 0j ) of the j *^ 
cap. 
(b) Patching together the Sj and Kj values for each spherical cap, one obtains two 
discrete functions S — S{6, (p) and K = K{6, (p) defined on the sphere S*^, which 
can be used as statistical indicators to measure the deviation from Gaussianity as 
a function of the angular coordinates {9, (f). The Mollweide projection of skewness 
and kurtosis functions S = S{9, (f>) and K = K{9, </>) are nothing but skewness 
and kurtosis maps (hereafter S*— map and /-iT— map). 

Clearly, the discrete functions S = S{9, (p) and K = K{9, (j)) can be expanded 
into their spherical harmonics in order to determine their power spectra Si and K^. 
Thus, for example, for the skewness one has S{9,(j)) = J27Lo^m=-£^^'m'^imid, (p) 
and Si = {2£ + 1)^^ J2m l^^ml^- Similar expressions obviously hold for the kurtosis 
K^Ki9,<P). 

In the remainder of this work we shall report the results of our Gaussianity 
analysis performed with S = S{9, (p) and K = K{9, (p) indicators calculated from 
single frequency and foreground reduced maps with and without a KQ75 mask. To 
minimize the statistical noise, in the calculations of 5*— map and iiT— map from the 
input maps, we have scanned the celestial sphere with spherical caps of aperture 
7 = 90°, centered at 12 288 points homogeneously distributed on the two-sphere. 

For the sake of brevity we do note show examples of S and K maps, which 
provide only qualitative information on large-angle deviation from Gaussianity (see 
figures in Refs. [SHH for examples). To obtain quantitative information about the 
large angular scale (low i) distributions for the non-Gaussianity S and K maps ob- 
tained from the CMB input maps used, we have calculated the (low £) power spectra 
Sg and Kg for S and K maps. The deviation from Gaussianity and the statistical 
significance were estimated by comparing these power spectra with the correspond- 
ing averaged power spectra Se and Ke calculated from S and K maps obtained 
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by averaging over 1 000 Monte-Carlo-generated statistically Gaussian CMB maps. 
To have an overall assessment of low i power spectra Sg and Kg calculated from 
each CMB input, we have performed a y^ test to find out the goodness of fit for 
Si and Kg multipole values as compared to the expected multipole values from the 
MC Gaussian maps. In this way, we obtained one number for each map that collec- 
tively ('globally') quantifies the deviation from Gaussianity. For the power spectra 
Sg calculated from 5*— maps obtained from the five- years maps with a KQ75 mask 
we found that the ratio x^/ dof (dof stands for degree of freedom) from the K, Ka, 
Q, V, and W maps are given, respectively, by 21.5, 4.9, 6.0, 5.2, and 3.9, while for 
the kurtosis power spectra Kg of these maps the values of x^/dof are, respectively, 
35 652, 135, 0.5, 6.4, and 5.6. Clearly a good fit occurs when x^/dof '^ 1. More- 
over, the greater are the x^/dof values the smaller are the x^ probabilities, that 
is the probability that the multipole values Sg and Kg and the expected MC mul- 
tipole values agree. For Sg and Kg obtained from the full-sky foreground-reduced 
five years ILC input maps we found that x^/dof are, respectively, 35.7 and 2 368. 
These results reduce to 1.2 and 0.4 when the KQ75 mask is employed. In brief, our 
analyses show that the unmasked band maps are significantly non-Gaussian but 
the deviation from Gaussianity is substantially reduced to a level compatible with 
Gaussianity for Q and V maps, whereas the full-sky foreground reduced five years 
ILC mask is again significantly non-Gaussian but the level of non-Gaussianity drops 
to a level that is consistent with Gaussianity when the KQ75 mask is used. 
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